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Abstract. A smooth curve of degree d in the real projective plane is hyper- 
bolic if its ovals are maximally nested, i.e. its real points contain [f J nested 
ovals. By the Helton- Vinnikov Theorem, any such curve admits a definite sym- 
metric determinantal representation. We use polynomial homotopy continua- 
tion to compute such representations numerically. Our method works by lifting 
paths from the space of hyperbolic polynomials to a branched cover in the space 
of pairs of symmetric matrices. 



Introduction 

Let p G C[t,x,y] be a homogeneous polynomial of degree d ^ 1. A (linear 
symmetric) determinantal representation of p is an expression 

p = det(tMi + xM 2 + yM 3 ) 

where Mi, M 2 , M3 are complex symmetric matrices of size d x d. Determinantal 
representations of plane curves are a classical topic of algebraic geometry. Exis- 
tence in the general case, at least for smooth curves, was first proved by Dixon in 
1902 [2]. For an exposition in modern language, see Beauville [I]. 

Real determinantal representations of real curves have only been studied sys- 
tematically much later in the work of Dubrovin |3] and Vinnikov |12| . Of particular 
interest here are the definite representations, where some linear combination of 
the matrices Mi, M 2 , M 3 is positive definite. By a celebrated result due to Helton 
and Vinnikov [6], these correspond exactly to the hyperbolic curves, whose real 
points consist of maximally nested ovals in the real projective plane. 

The Helton- Vinnikov theorem (previously known as the Lax Conjecture) has 
attracted attention in connection with semidefinite programming, since it charac- 
terizes the convex subsets of the real plane that can be described by linear matrix 
inequalities. See Vinnikov |14j for an excellent survey. 

While the Helton- Vinnikov theorem ensures the existence of a definite determi- 
nantal representation for any hyperbolic curve, finding such a representation for 
a given polynomial p remains a difficult computational problem. With a suitable 
choice of coordinates, we can restrict to representations of the form 

p = det(tl d + xD + yR) 

where is the identity matrix, D is a real diagonal and R a real symmetric matrix. 
The hyperbolicity of p is reflected in the fact that for any point (u, v) G M 2 , all 
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roots of the univariate polynomial p(t, u, v) G M[t] are real. Given such p, the 
computational task of finding the unknown entries of D and R leads, in general, 
to a zero-dimensional system of polynomial equations. However, as d grows, this 
direct approach quickly becomes infeasible in practice. This, as well as symbolic 
methods and an alternative approach via theta functions based on the proof of 
the Helton- Vinnikov theorem, have been investigated in (TUJ. As far as actual 
computations are concerned, d = 6 was the largest degree for which computations 
terminated in reasonable time. 

In this note, we present a more sophisticated numerical approach, implemented 
using the NumericalAlgebraicGeometry package [7] in MACAULAY2 [I] (NAG4M2). 
We consider the branched cover of the space of homogeneous polynomials by pairs 
of matrices (D, R) (with D diagonal and R symmetric) via the determinantal map 
(D, R) i — y det(tId + xD + yR). We use known results on the number of equivalence 
classes of complex determinantal representations to show that the determinantal 
map is unramified over the set of smooth hyperbolic polynomials (Thm. II. 3p . We 
then use the fact that this set is path-connected. In fact, an explicit path connect- 
ing any hyperbolic polynomial to a certain fixed polynomial was constructed by 
Nuij in P|, which we refer to as the N-path. Our algorithm works by constructing 
a lifting of the N-path to the covering space. The advantage over just solving the 
zero- dimensional system of equations directly via homotopy continuation is that 
we need only track a single path, instead of one path for each solution. We also 
recover the Helton- Vinnikov theorem from these topological considerations and 
the count of equivalence classes of complex representations. 

At this point, we do not know whether the N-path indeed avoids the singular 
locus. We prove this for the case d = 2 by considering the corresponding discrim- 
inant (Prop. 12. 2h . At any rate, since the singular locus has codimension at least 
2 inside the set of strictly hyperbolic polynomials, the N-path can be expected to 
avoid singularities for almost all of the starting hyperbolic polynomials. 

To account for the unlikely scenario that the (real) N-path goes through the 
singular locus, we provide an algorithm that produces a complex determinantal 
representation via the complexification of the N-path. The fact that this repre- 
sentation is real with a non-zero probability gives a straightforward probabilistic 
algorithm for finding a real determinantal representation. 

Our current proof-of-concept implementation uses standard double floating 
point precision and evaluation of polynomials in the top-level interpreted language 
of MACAULAY2. Even with these limitations we can compute some examples for 
d = 7 within one day. With arbitrary precision arithmetic and speeding up the 
numerical evaluation procedure, we see no obstacles to computing robustly for d 
in double digits using the present-day hardware. 

Acknowledgements. We are grateful to Institut Mittag-Leffler, where this project 
has started, for hosting us in the Spring of 2011. We would like to thank Greg 
Blekherman for helpful discussions concerning the discriminant computed in Sec- 
tion [2 
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1. DETERMINANTAL MAPPING 



For d ^ 1, consider the affine space of polynomials 



T = {p G C[£, x, y] \p is homogeneous of total degree d and p(l, 0, 0) = l} 



and the vector space of pairs of complex resp. real symmetric matrices 



M = {(D,R)e (Sym d (C)) 2 | D is diagonal} 
M u = M n (Sym(R)) 2 . 



The dimensions of these four spaces (over C resp. R) are all the same, namely 



and its restriction to .Mr. A polynomial p G J-r is called hyperbolic (with respect 
to (1,0,0)) if all roots of the univariate polynomial p(t,u,v) G R[t] are real, for 
all (u,v) G R 2 . It is called strictly hyperbolic if the roots of p(t,u,v) are real and 
distinct, for all (u, v) G R 2 . Write 



The image of Mr under $ is contained in 7i. It is also not hard to show that it 
is closed (see [TT| Lemma 3.4]). Our first goal is to find a connected open subset 
U of H such that the restriction of $ to <3> _1 (£7) is smooth. 

Proposition 1.1. 

(1) The set % is closed in J-r ; contractible and path-connected. 

(2) The interior int('H) ofH is the set of strictly hyperbolic polynomials and 
% is the closure ofmt(H). 

(3) A polynomial f G Ti is strictly hyperbolic if and only if the projective 
variety Vc(f) defined by f has no real singular points. 

(4) Let H° be the set of hyperbolic polynomials f G H for which Vc(jf) is 
smooth. Then int("H) \ W has codimension at least 2 in J-r. 

Proof. (1) and (2) are proved by Nuij [U] (see also below). (3) is proved in [TTj 
Lemma 2.4]. (4) follows from the fact that the elements of int('H) have no real 
singularities, while complex singularities must come in conjugate pairs. □ 

For fixed p G %, the group SL^(C) acts on the determinantal representations 
p = det(tMi + xM 2 + yM 3 ) via symmetric equivalence. In other words, any 
A G SL d (C) gives a new representation p = det(L4MiA T + xAM 2 A T + yAM 3 A T ). 
When we restrict to the normalized representations we are considering, we have 
an action on pairs (D,R) G $ _1 (p) by those elements A G SL^(C) for which 
AA T = I d (i.e. A G SO d (C)) and ADA T is diagonal. 



We study the map 




M -> T 

(D, R) i-> det(tI d + xD + yR). 



H = {p G Jk | p is 



hyperbolic}. 
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Theorem 1.2. Any p G J 7 has only finitely many complex representations p = 
det(tld + yD + zR) up to symmetric equivalence. If the curve Vc{p) is smooth, 
the number of equivalence classes is precisely 2 9 ~ 1 ■ {2 9 + 1), where g = ( ~ x ) is 
the genus ofVcip)- 

Proof. For smooth curves, the equivalence classes of symmetric determinant al rep- 
resentations are in canonical bijection with ineffective even theta characteristics; 
see |1U1 Thm. 2.1] and references given there. □ 

Theorem 1.3. The set W of smooth hyperbolic polynomials is an open, dense, 
path- connected subset of 7-L, and each fibre of $ over a point of 7i° consists of 
exactly 2 9_1 ■ (2 9 + 1) • 2 d_1 • d\ distinct points. 

Proof. Let p G int("H) and let (D, R) G $ _1 (p), which means p = det(tld + xD + 
yR). The diagonal entries of D are the zeros of p(t, — 1, 0). Since p is strictly 
hyperbolic, these zeros are real and distinct. So D is a real diagonal matrix 
with distinct entries. It follows then that the centralizer of D in SOd(C) consists 
precisely of the 2 d diagonal matrices with entries ±1. Let S be such a matrix with 
S 7^ ±1. We want to identify the set of symmetric matrices R that commute with 
S. Up to permutation, we may assume that the first k diagonal entries of S are 
equal to —1 and the remaining d — k are equal to 1. It follows then that any R with 
SR = RS must have = rji — if % > k ^ j, so that R is block-diagonal. For 
such R to show up in a pair (D, R) G the polynomial p must be reducible. 

In particular, if p G "H°, there is no such S commuting with R. It follows then 
that {SRS | S diagonal with S 2 = 1} has 2 d ~ l distinct elements. Permuting the 
distinct diagonal entries of D gives d\ possible choices of D. This, combined with 
the count of equivalence classes in the preceding theorem, completes the proof. □ 

Corollary 1.4. The restriction o/$ to $~ 1 ('H ) is smooth. 

Proof. The restriction of $ to $~ 1 ('H°) is a polynomial map with finite fibres that 
is unramified over % , since the cardinality of the fibre does not change. Hence it 
is smooth (see for example Hartshorne [HI III. 10]). □ 

Corollary 1.5 (Helton- Vinnikov Theorem). Every hyperbolic curve possesses a 
definite real determinantal representation. 

Proof. Since all fibres of $ over H° have the same cardinality and "H° is path- 
connected, the number of real points in each fibre must also be constant over W . 
That number cannot be zero, since there clearly exist fibres with real points. It 
follows that W is contained in $(.Mir). On the other hand, $(A^ir) is closed in 
J-jj by [TT] Lemma 3.4] and contained in H, hence $(A^u) = H. □ 

Remark 1.6. The number of equivalence classes of real definite representations 
of a hyperbolic curve is in fact also known, namely it is 2 9 . See [TU] and references 
to [13] given there. We conclude that $ _1 (p) D .Mr consists of 2 9 ■ 2 d ~ 1 ■ d\ distinct 
points for every p G TL°. 

Note also that, even if p is hyperbolic, it will typically admit real determinantal 
representations det(tM! +yM 2 + zM 3 ) that are not definite, i.e. are not equivalent 
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to such a representation with Mi = Id and M 2 , M 3 real. Such representations do 
not reflect the hyperbolicity of p. 



2. NUIJ-PATH 

In order to use homotopy continuation methods for numerical computations, we 
need an explicit path connecting any two given points in the space % of hyperbolic 
polynomials. Following Nuij [9], we consider the operators 

dp 



T k , s ■ P >-> p + sx k 
G s : p\-> p(t, sx) 



Of 



rp rpd rpd 

r s — 1 l jS ' ' ' 1 n,s- 

Proposition 2.1 (Nuij [9]). Forp G mt(/H), we have H s (p) = Fi_ s G s (p) e int('H) 
for all s G [0, 1], with FqGi(p) = p and FiGo(p) not depending on p. □ 

We call [0, 1] 3 s i— > H s (p) the N-path of p. In order to ensure smoothness of 
the map $: Ai — > % along an N-path, we would like to ensure that the N-path 
stays inside the set "H of smooth hyperbolic polynomials and thus away from the 
ramification locus. 



Proposition 2.2. Let d 

se [0,1). 



2. Then H s {p) G W holds for all p G H and all 



Proof. Let D = Diag(di, d%) and R 



. The quadric 



r (i,i) r (i,2) 

7(1,2) r( 2 ,2)_ 

H S {${D, R)) = Ax 2 + Bxy + Cy 2 + Dxt + Eyt + Ft 2 
is not contained in T-L° if and only if it factors, which happens if and only if 

" A B/2 D/2 



4Disc(H s (<S>(D,R))) = 4det 



B/2 C E/2 
D/2 E/2 F 



2s 2 {s 
2s 2 (s 

,4„2 



2 ^ d 2 ) 2 



D 2 {di 



I) 2 (m 



r 2 2 



s % (di - d 2 y + 



8r 2 12 s 2 (1 



,) 2 + 



16(s- l) 4 
= 0. 

The sum-of-squares representation was produced using the intuition obtained by 
a numerical sum-of-squares decomposition delivered by YALMIP and further ex- 
act symbolic computations in MATHEMATICA (see Nuij -d2 .mathematica [8]). 
The components of the sum-of-squares decomposition above vanish simultaneously 
only when s — 1 (this proves the proposition) and either rvi — or d 2 — di. □ 
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Remark 2.3. There is numerical evidence that in case d = 3 the N-path stays 
away from the locus singular cubics. 

An attempt to plug in H S (§(D, R)) into the discriminant freezes all computer 
algebra systems: the result, a polynomial of degree 36 in 9 variables and di, is 
an extremely long expression if written in the monomial basis. 

While it seems infeasible to derive symbolically that (T-L\H°) D {N-path} has no 
real points except for s = 1 in a way similar to Proposition 12. 2} our experiments 
driven by numerical real optimization in MATHEMATICA support this conclusion 
(see Nuij -d3 .mathematica [S]). 

In view of the preceding proposition, it seems natural to ask the following. 

Question 2.4. It it true for all p that the N-path H s (p), with s G [0, 1), does not 
meet the ramification locus of the cover map $? 

Should the answer be no, one can still attempt to avoid the ramification points in 
the following manner. 

Proposition 2.5. The N-path H s (p), with parameter s varied along a piecewise- 
linear path [0, c] U [c, 1) C C, does not meet the ramification locus of $, for almost 
all c & C (More precisely, this holds for any c taken in the complement of some 
proper real algebraic subset of C ~ M 2 .) 

Proof. This immediately follows from the fact that the ramification locus is a 
proper complex subvariety of J 7 and therefore has real codimension at least 2. □ 

Remark 2.6. Note that using a random path as described in the Proposition may 
result in a non-real determinantal representation: indeed, a non-real path for s is 
not guaranteed to result in a real point in the fiber $~ 1 (iJ (j*)) when a real point in 
is taken. Nevertheless, the probabilistic algorithm suggested by the 
Proposition gives a way to construct not only complex, but real determinantal 
representations. 

One can ask for the probability of obtaining a real determinantal representation 
at the end of the path described in Proposition 12.51 (For a more precise question, 
one may pick c on a unit circle with a uniform distribution.) While this probability 
is clearly non-zero, deriving an explicit lower bound seems to be a very hard 
problem. A naive intuition suggests that the probability can be estimated as a 
ratio of the number of real representations (Remark ll.6p to the total count of 
complex representations (Theorem II. 31) . The experiments with quartic hyperbolic 
curves, however, suggest that the probability is much higher. 

3. Algorithm and implementation 

Given a hyperbolic polynomial p G H, the N-path H s (p) connects p = Hi{p) 
with Pq = H (p) which does not depend on p. This suggests the following algo- 
rithm to compute a determinantal representation (D p , R p ) G .Mr for p: 

(1) Pick (D q , R q ) G Aim. giving a strictly hyperbolic polynomial q = &{D q , R q ). 
Track the homotopy path H s (q) from s = 1 with the start solution (D q , R q ) 
to s = producing the target solution (D po , R Po ). Then p = $(_D po , R po ). 
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(2) Track the homotopy path H s (p) from s = with the start solution (D po , R po ) 
to s — 1 to obtain (D p , R p ) such that p = $(-D p , R p ). 

Note that, in principle, the first step only has to be performed once in each 
degree d. In what follows we describe two ways to set up a polynomial homotopy 
continuation for the pullback of an N-path. 

3.1. N-path in the monomial basis. One way is to take the coefficients of the 
polynomial <&(D,R) — H s (p) G C[D, R,s][x,y,t] with respect to the monomial 
basis of J 7 . This gives a family of square (^equations= ^unknowns) systems of 
polynomial equations in C[D,R] parametrized by s. 

Then this family is passed to a homotopy continuation software package (we 
use NAG4M2 [7]). As long as s G C follows a path that ensures that H s (p) stays 
in "H° C H, there are no singularities on the homotopy path except, perhaps, at 
the target system (see discussion in Section [2D- 

The bottleneck of this approach is the expansion of the determinant in the 
expression $>(D,R) and evaluation of its (x, y, ^-coefficients: it takes Q(d\) oper- 
ations and results in the expression with Q(d\) terms. This limits us to d < 5 in 
the current implementation of this approach. 

3.2. N-path with respect to a dual basis. While it may seem that picking a 
basis of J 7 different from the monomial one does not bring any advantage, it turns 
out to be crucial for practical computation in case of larger d. 

We fix a dual basis in J 7 * consisting of m — dim J 7 evaluations e, at general 
points (xi,yi,ti) G C 3 , for i = l,...,m. The current implementation generates 
the points with coordinates on the unit circle in C at random. 

Now the family of polynomial systems to consider is 

fH = e^D, R) - H s {p)) G C[D, R, s], i = 1, . . . , m. 

Since ei(<&(D, R)) = det(Jtj + Dxi + Ryi), the evaluation of hi and its partial 
derivatives costs 0(d 3 ). We modified the NAG4M2 implementation of evaluation 
circuits (a.k.a., straight-line programs) to include taking a determinant as a basic 
evaluation node. 

4. Example 

The last improvement in the implementation allows us to compute examples for 
larger d. With an implementation of the homotopy tracking in arbitrary precision 
arithmetic, we see no obstacles to computing determinantal representations for d 
in double digits. 

To give an example, we choose the sextic 

p = - 36x 6 - 157xy - 20xV - 109x 2 y 4 + 246xy 5 - 92y 6 - 12x 3 y 2 t + 90x 2 y 3 t 
+ lOxyH + 76yH + 49xH 2 + 156x 2 y 2 t 2 - 16xy 3 t 2 + 132yH 2 + 12xy 2 t 3 
- 14y 3 t 3 - 14x 2 t 4 - 27y 2 t A + t 6 . 



8 



ANTON LEYKIN AND DANIEL PLAUMANN 



The polynomial p is hyperbolic, since p = $(-D, R) with 
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(4.1) D = Diag(— 3, —2, —1, 1, 2, 3), R 



Assuming this pair (D, R) is not known, let us describe the application of our 
algorithm to recover a determinantal presentation of p; one can reproduce the fol- 
lowing results by running lines in showcase .m2 [S]. First, taking an arbitrary pair 
(D q ,R q ) and tracking the N-path H s (q) from the strictly hyperbolic polynomial 
q = $(-D g , R q ) = Hi(q) to the fixed polynomial po = Ho(q) we get 

D P0 = Diag(.222847, 1.18893, 2.99274, 5.77514, 9.83747, 15.9829) 



R 



6 

2.51352 
1.19571 
4.04309 
1.42786 
-1.98597 



2.51352 
6 

3.08656 
.468873 
2.38468 
1.05948 



1.19571 
3.08656 
6 

.785785 
4.66027 
2.29433 



4.04309 
.468873 
.785785 
6 

1.6226 
.933245 



1.42786 
2.38468 
4.66027 
1.6226 
6 

3.50198 



-1.98597' 
1.05948 
2.29433 
.933245 
3.50198 
6 



Hq(p) = H (q) to p = Hi(p), we obtain 



Tracking the N-path H s (p) from p 

£>' = Diag(-3, -2, -1,1,2,3) 

.596508 -1.43241 2.00316 1.10471 -.725394' 

.596508 .739773 1.79407 .0604427 -1.60948 

-1.43241 .739773 1.56816 1.66137 -.165953 

2.00316 1.79407 1.56816 .839374 2.00885 

1.10471 .0604427 1.66137 .839374 



R' 



-.725394 -1.60948 



1.57679 
.165953 2.00885 1.57679 



which is an alternative determinantal presentation of p. While we returned to the 
same point p in the base of the cover $, the route taken has led us to a different 
sheet than the sheet of the fiber point (D, R) e $ _1 (p) in (14. ip used to construct 
this example. 

With the default settings of NAG4M2 the homotopy tracking algorithm takes 
28 steps on the first path and 15 steps on the second. 
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